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Abstract — In  this  paper  we  present  a  strategy  that  allows 
a  swarm  of  autonomous  tugboats  to  cooperatively  move  a 
large  object  on  the  water.  The  two  main  challenges  are:  (1) 
the  actuators  are  unidirectional  and  experience  saturation;  (2) 
the  hydrodynamics  of  the  system  are  difficult  to  characterize. 
The  primary  theoretical  contribution  of  the  paper  addresses 
the  first  challenge.  We  present  a  tracking  controller  and  force 
allocation  strategy  that,  despite  actuator  limitations,  result 
in  asymptotically  convergent  tracking  for  a  certain  class  of 
reference  trajectories.  The  primary  practical  contribution  is 
the  introduction  of  a  set  of  adaptive  control  laws  that  address 
the  second  challenge  by  compensating  for  unknown,  and 
difficult  to  measure,  hydrodynamic  parameters.  Experimental 
verification  of  the  controllers  is  presented  using  a  1:36  scale 
model  of  a  U.S.  Navy  ship,  inside  the  Naval  Academy’s  unique 
380  ft  testing  tank. 

I.  Introduction 

ROBOT  swarms  are  large  groups  of  small,  rela¬ 
tively  unsophisticated,  robots  working  in  concert 
to  achieve  objectives  that  are  beyond  the  capability  of 
a  single  robot.  One  example  of  an  application  that  can 
benefit  from  this  approach  is  non-prehensile  cooperative 
manipulation,  where  a  group  of  non-articulated  mobile 
robots  attempts  to  transport  a  larger  object  in  the  plane, 
by  applying  forces  to  its  perimeter.  The  advantages 
of  the  swarm  are:  (1)  its  ability  to  distribute  applied 
forces  over  a  large  area,  achieving  an  enveloping  grasp 
on  large  objects;  and  (2)  the  maximum  wrench  the 
swarm  can  exert  increases  linearly  as  the  number  of 
swarm  members  increases.  We  are  particularly  interested 
in  marine  applications  involving  autonomous  tugboats 
such  as  towing  disabled  ships  (ex.  U.S.S.  Cole),  trans¬ 
porting  components  of  large  offshore  structures  (ex.  oil 
platforms),  or  positioning  littoral  protection  equipment 
(ex.  hydrophone  arrays).  However,  most  of  our  group’s 
work  can  be  extended  to  ground  robots  as  well. 

Swarm  control  strategies  represent  a  very  active  area 
of  research  (see  for  example  [14],  [25],  [7],  [4]  and 
references  within).  However,  most  of  those  works  only 
consider  kinematic  objectives  (i.e.,  position  and  velocity), 
such  as  collective  motion  (flocking)  or  mapping  tasks. 
In  cooperative  manipulation,  second  order  forces  and 
contact  mechanics  must  be  taken  into  account.  Such 


issues  have  been  examined  before  in  the  context  of  single 
robots  or  small  groups  of  two  to  three  robots  [12],  [18], 
[24],  [22],  Behavior  based  approaches  yield  interesting 
results  but  lack  performance  guarantees  [20],  [15],  [8], 
[H]- 

Marine  applications  with  a  high  number  of  actuators 
have  been  considered  before  as  well,  but  exclusively  in 
the  context  of  large  vessels  with  “thruster  pods”  such  as 
commercial  cruise  ships  [6],  [9],  [17]  and  [27].  As  such, 
in  each  of  these  works,  the  control  scheme  is  completely 
centralized  with  no  uncertainty  in  the  location  or  actions 
of  the  actuators. 

Closely  related  works  on  distributed  swarm  manip¬ 
ulation  are  [23],  [26],  and  [16].  Within  these  works, 
controllers  were  designed  to  force  robots  to  surround 
the  object.  Inter-robot  spacing  was  constrained  to  be 
small  enough  that  it  is  impossible  for  the  object  to 
“escape”  meaning  that  as  the  robots  move,  so  must  the 
object.  While  this  approach  is  decentralized,  the  primary 
drawback  is  that  it  was  designed  for  land-based  robots. 
Therefore,  they  treat  the  task  as  a  position  control 
problem  and  ignoring  the  dynamic  forces,  momentum 
and  actuator  saturation. 

In  this  paper  we  present  a  control  strategy  that  enables 
the  barge  to  track  a  reference  trajectory.  It  accounts  for 
second  order  dynamics,  actuator  limitations,  and  model 
uncertainty  associated  with  linear  drag  effects.  A  special 
property  of  the  controller,  when  used  in  conjunction 
with  the  force  allocation  strategy  presented  here  and 
grasp  synthesis  approach  in  companion  work,  is  that 
it  is  asymptotically  stable  despite  actuator  limitations 
on  individual  swarm  members.  In  the  authors’  previous 
work  [2],  a  commutation  strategy  was  employed  to 
exactly  distribute  the  force/torque  command  among  the 
swarm  members;  however,  the  tugboats  must  be  aligned 
in  a  precise  opposing  pair  configuration  along  the  barge 
hull.  Though  this  approach  offers  a  method  by  which  to 
guarantee  that  the  swarm’s  combined  effort  produces  the 
commanded  request;  maintaining  the  exact  alignment 
may  be  difficult  in  heavy  sea  states.  Therefore  in  this 
work,  we  wish  to  investigate  the  possibility  of  utilizing 
an  alternative  optimization  approach  for  force/torque 
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allocation.  In  addition,  in-water  experiments  show  that 
it  is  also  able  to  compensate  for  model  uncertainties 
associated  with  drag  coefficient  uncertainties,  which  is  a 
critical  issue  for  practical  implementation  in  the  marine 
setting. 

The  remainder  of  the  paper  is  organized  as  follows. 
In  Section  II,  we  provide  the  second  order  dynamic 
model  of  the  marine  vessel  and  formalize  the  problem 
statement.  In  Section  III,  we  present  the  high  level  track¬ 
ing,  adaptive  controller  along  with  our  force  allocation 
methodology  in  Section  IV.  In  Section  V,  simulation 
results  are  provided,  and  the  experimental  test  bed  and 
results  are  presented  in  Sect.  VI.  Section  VII  summarizes 
our  contributions  and  discusses  our  future  work. 


II.  System  Model 

We  refer  to  the  object  to  be  manipulated,  as  a  barge. 
Mathematically,  it  is  defined  by  a  compact,  convex  set 
O  C  R2  .  Define  a  body-fixed  coordinate  frame  (i,  j) 
attached  to  the  centroid  of  O  (see  Figure  1).  Ignoring 
heave  (vertical),  pitch  and  roll  motions,  a  planar  3  DOF 
model,  low  speed  model  is  [6] 

P  =  R(ip)v  (1) 

Mv  +  Dv  =  Unet  (2) 


where  P  =  [x,y,ip]T  denotes  the  inertial  ( x,y )  position 
and  barge  heading  (ip),  v  =  [vx,vv,ix>]T  is  the  barge’s 
velocity  expressed  in  the  body  frame  (also  called  surge, 
sway  and  yaw  rate).  R(ip)  denotes  the  rotation  matrix 
from  the  body  frame  to  the  inertial  frame  and  is  defined 
as  follows 


R(iP) 


cos  (ip)  —  sin(V’)  0 
sin(V’)  cos  (ip)  0 

0  0  1 


(3) 


M  =  diag {mx,my,  Iz}  €  R3x3  represents  the 
mass/inertia  matrix  that  includes  added  mass  [6],  and 
D  =  diag {XU,YV,  Nr}  €  R3x3  represents  the  hydrody¬ 
namic  drag  matrix  (for  our  low  speed  application,  only 
linear  drag  is  considered). 

The  right  hand  side  of  (2)  is  the  net  applied  wrench  on 
the  barge,  unet  =  [Fx,Fy,Mz]T  due  to  the  N  tugboats. 
Note  that  in  practice,  actual  tugboats: 

1)  “tie-up”  to  the  barge,  meaning  that  their  configu¬ 
ration  is  nominally  time  invariant  and  their  ability 
to  transfer  forces  is  not  friction  dependant; 

2)  have  propellers  and  engines  that  are  significantly 
less  efficient  when  run  in  reverse;  and 

3)  have  upper  limits  on  the  thrust  they  can  apply 
very  often  two  or  more  tugs  are  required  to  generate 
enough  force  to  move  the  barge. 

Reflecting  these  assumptions,  unet  =  Bu,  where 
u  =  [ui, . .  .un]t  G  U  is  a  vector  of  input  push  force 
magnitudes  for  each  tugboat.  The  set  U  is  described 
as  0  <  Ui  <  wmax,V?'  G  [1,TV].  B  G  R3xJV  captures 
the  geometric  configuration  of  the  tugs  and  requires 
additional  notation. 


Fig.  1.  Swarm  Manipulation  Scenario.  N  robots(dark  circles) 
are  attached  to  the  object  (shaded  polygon);  and  can  apply  forces 
at  some  fixed  incident  angles.  The  position  of  the  ith  boat  is 
determined  by  angle  6i  and  the  push  direction  by  the  angle  oli  . 


The  barge  boundary  is  denoted  by  dO ,  assuming  it 
is  convex  allows  us  to  parameterize  the  points  on  the 
boundary  using  the  angle  9 ,  measured  counter  clockwise, 
relative  to  the  x  axis  of  the  body  frame: 

{x(9)i  +  y(9)j  G  dO,  ^9  G  S1}  (4) 


in  addition,  define  a  parameter  a  G  S1  which  indicates 
the  direction  of  the  push  force,  measured  counter  clock¬ 
wise  relative  to  the  i-axis  of  the  body  frame.  Then  the 
contact  configuration  of  the  ith  robot  can  be  described 
by  a  vector  [0j,cq]T .  The  ith  column  is 


Bi  = 


COs(cq) 

sin(cq) 

— y(9i )  cos  cq  +  x(9i)  sin  cq 


(5) 


Note  that  there  is  no  constraint  on  the  allowed  push 
directions,  cq,  as  there  would  be  in  the  case  of  friction 
assisted  grasping.  Throughout  this  paper  we  assume  that 
the  configuration  of  the  tugboats  (and  hence  the  B 
matrix)  represents  a  force  closure  grasp.  Force  closure 
is  defined  as  the  ability  of  a  robot  grasp  to  resist  or 
apply  a  wrench  in  an  arbitrary  direction.  Due  to  the 
unilateral  constraint  that  the  tugboat  forces  u  are  non¬ 
negative,  in  order  to  achieve  force  closure  the  set  of 
possible  applied  wrenches  at  all  the  contacts  (i.e.  the 
columns  of  B)  must  positively  span  R3;  this  condition  is 
easily  tested  [13].  For  this  application  stating  that  the 
grasp  is  a  force  closure  configuration,  is  equivalent  to 
saying  that  the  position  and  orientation  of  the  object 
is  small  time  locally  controllable,  under  the  constraint 
that  the  applied  forces  are  nonnegative.  It  is  a  well  known 
result  that  at  least  4  contacts  [13]  are  required  to  satisfy 
this  condition  for  the  class  of  objects  and  contact  types 
considered  here.  Therefore,  in  this  paper  we  assume  that 
the  swarm  size  TV  >  4.  Note  that,  given  the  B  matrix  it 
is  possible  to  compute  an  upper  bound  on  unet  consistent 
with  the  constraints  on  u 


'Uj\ net 


min 

jkeNxN 


^max 


N 

max(0, 

2  —  1 


Bj  x  Bk 
II Bj  x  Bk || 


•  Bi). 


(6) 
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Synthesizing  the  optimal  placement  of  the  tugs  is  beyond 
the  scope  of  this  paper,  but  is  addressed  in  other  work 
from  our  group  [3]. 

III.  Tracking  Control  Design 
Problem  3.1:  Given,  the  nominal  ship  model  from  (1), 
a  reference  trajectory  Pd(t),  measurements  of  P  and 
P,  and  the  upper  bound  on  the  thrust  the  swarm  can 
generate  in  all  directions  unet\  determine  a  feedback  law 
U  :  R6  ->  unet  such  that  ||unet(t)||  <  unet,Vt  £  R+. 

In  order  to  promote  this  position  objective,  the  fol¬ 
lowing  error  signals  are  defined 

ep  =  PT0 <p)  ( Pd~P ) ,  ev  —  v  —  RT(i’)Pd  (7) 

where  Pd(t)  £  R3  represents  a  sufficiently  smooth  (i.e. , 
differentiable)  position/orientation  trajectory.  In  order 
to  facilitate  the  design  of  the  2nd  order  system,  the  filter 
tracking  error  r(t)  £  R3  is  specified  as  follows 


r(t)  =  e„  +  kpep  (8) 

where  kp  £  R+  denotes  a  positive,  scalar  control  gain. 
In  order  to  account  for  unknown  hydrodynamic  drag 
parameters,  the  following  parameter  estimation  error 
signal  0(f)  £  R3  is  defined  in  the  following  manner 

0  =  0-0  (9) 

.  iT 


where  0(f)  = 


Xn,  Y„,  Nr 


denotes  the  yet  to  be 


designed  parameter  estimate  vector.  In  order  to  proceed 
with  the  design  of  the  tracking  controller,  we  first  develop 
the  open-loop  position/ velocity  tracking  systems  [19] 


ep  =  -S(v)ep  +  e„ 

e„  =  v  +  S(v)RT(iP)Pd-RT(ip)Pd 

where  the  skew-symmetric  matrix  S(v)  is  given  as 

0  —  U!  0 

S(u)  =  |  w  0  0 

0  0  0 


(10) 


(11) 


In  order  to  design  the  control  input  unet,  the  open-loop 
filtered  tracking  error  dynamics  are  developed  in  the 
following  manner 

Mr  =  [-Y  (i/)0  +  unet]  +  M  S(u)  (. RT(ip)Pd  -  kpe^j 
-RT(ip)Pd  +  kpev 

(!2) 

where  the  regression  matrix  Y  (v)  and  the  unknown  pa¬ 
rameter  vector  0  are  defined  by  the  following  expressions 


1»  = 


u 

0 

0  ' 

'  xu  ' 

0 

V 

0 

,0  = 

Yv 

0 

0 

LO 

Nr 

(13) 


Based  on  the  open-loop  tracking  error  dynamics,  the 
vessel  control  input  unet(t)  to  promote  position  tracking 
is  specified  in  the  following  manner 


'Unet 


=  Y(v)Q  —  M  S(v)  [rt  (ip)Pd  -  kpe^ 


-RT  (ip)  Pd  +  kpev 


-ep-  krr 


(14) 


where  kr  £  R+  denotes  a  positive,  scalar  control  gain 
(note  that  kr  could  be  defined  as  a  positive  definite, 
diagonal  gain  matrix). In  order  to  calculate  a  suitable 
upper  bound  on  the  required  forces/torque  of  (14),  the 
following  projection  based  update  law  is  employed  to 
generate  a  bounded  parameter  estimate  vector  0(f) 

r  0  if  0  =  0,  —YT(v)r  >  0 

0  =  <  0  if  0  =  0,  —YT(v)r  <  0  (15) 

[  —YT(v)r  otherwise 

where  0  denotes  the  upperbound  values  for  the  drag 
coefficients  (assumed  known).  After  substituting  in  the 
control  force/torque  vector  unet(t)  into  the  open-loop 
dynamics  of  (12),  the  following  closed-loop  dynamics  are 
obtained  for  r(f) 

Mr  =  —Y(v)Q  —  ep~krr  (16) 

Stability  and  the  position  tracking  result  can  be  il¬ 
lustrated  through  Lyapunov  analysis  by  defining  the 
following  non-negative  scalar  function  V  (f)  as 

V  =  \e?p  ep  +  |©T©  +  \rTMr  (17) 

After  taking  the  time  derivative  of  (17)  and  substituting 
in  the  expressions  of  (10),  (15),  and  (16),  the  time 
derivative  for  V  (f)  can  be  upperbounded  by  the  following 
expression 

V  <  -a||ep||2  -  Mr||2  (18) 

where  the  fact  that  e„  =  r  —  aep  has  been  utilized. 
Barbalat’s  Lemma  [10]  can  now  be  applied  to  (18)  to 
illustrate  asymptotic  position  tracking  in  the  sense  that 
lim^oo  ep(t),  r(t)  =  0. 

Remark  3.2:  In  an  effort  to  promote  a  more  fieldable 
solution,  one  of  the  objectives  of  this  work  was  to 
consider  actuator  constraints.  In  the  above  tracking 
controller,  the  commanded  force/torque  vector  unet(t ) 
is  not  directly  saturated  as  was  done  in  [5].  Rather,  the 
asymptotic  stability  result  coupled  with  an  appropriate 
desired  trajectory  generator  can  offer  some  ability  to 
influence  the  limits  on  unet(t).  Since  V(t)  <0  from  (18), 
then  it  is  known  that  after  some  time  the  error  signals 
may  be  neglected  from  the  unet(t)  leaving  the  following 
terms 

unet  »  Y(RTPd)e  +  MS(RTPd)RTPd  -  MRTPd  (19) 


where  we  have  utilized  the  assumption  that  enough 
time  has  elapsed  such  that  v  — >  RT (ip d)Pd-  From  the 
expansion  of  (19),  the  following  bounding  expression  can 
be  obtained 


l^net] 


XTl 


ma 


+  my\ipd\j  (\xd\ 

Nf\ipd\  +  Iz\i>d\ 


(l*d|  +  \Vd\)  ■ 
\Vd\)  H 


-mx(\xd\  +  \yd\) 
m.y(\xd\  +  \yd\) 


From  the  structure  of  the  above  expression,  the  selection 
of  the  desired  velocity/ accelerations  can  be  done  in 
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Fx  Force 


such  a  manner  as  to  ensure  that  the  command  vector 
Unet(t)  does  not  exceed  the  cooperative  capabilities  of 
the  swarm. 


IV.  Thrust  Allocation  Strategy 

With  the  design  of  unet(t ),  the  subsequent  task  is  to 
now  specify  the  individual  tugboat  thrusts  such  that 
their  combined  effort  generates  as  closely  as  possible 
unet  (t)-  Therefore,  the  thrust  allocation  strategy  must 
solve  the  following  linearly  constrained  least  squares 
problem. 

Problem  4.1:  Given  a  unet  from  (1),  determine  u  £ 
according  to 

min  \\ Bu  ~  unet\\2  (20) 

U 

such  that 


M  Torque 


0  - 

-1  ■ 


0  20  40  60  80  100  120 

Time  (sec) 


- Desired 

- -Optimization 

— —  Commutation 


- Desired 

-  -  -  Optimization 
- Commutation 


- Desired 

- -Optimization 

-  Commutation 


0  <  Ui  <  Mmax,  i  =  l,...,N.  (21) 

Remark  4.2:  The  linearly  constrained  least  squares 
problem  is  known  to  be  convex,  therefore  it  returns  a 
globally  optimal  solution.  Furthermore,  its  solution  is 
efficient  enough  to  be  computed  in  real  time  as  part  of 
the  feedback  law  and  scales  well,  0(N2),  with  swarm 
size  [1]. 


V.  Simulation  Results 


The  primary  focus  of  the  simulation  studies  is  to 
investigate  the  performance  of  the  optimization  approach 
of  Section  IV  for  force/torque  allocation.  To  this  end,  a 
barge  having  the  following  mass  and  drag  values  [17]  was 
simulated  within  the  MATLAB  SIMULINK  environment 

M  =  diag{  19.0  (kg),  35.2  (kg),  4.2  (kg  ■  m2)  } 

D  =  diag{  4.0  (££*),  10.0  (^),  1.0  (Nm  ■  sec)  } 

(22) 

The  swarm  configuration  about  the  barge  hull  was 
selected  according  to  [19]  as  follows 


ai  =  180°, 

a2  =  270°, 
a3  =  270°, 
aq  =  0°, 
a5  =  90°, 
a6  =  90°, 


rq  =  0.60  (m), 
r2  =  0.27  (m), 
r3  =  0.19  (m), 
r4  =  ri, 
rb  =  r2, 

=  r3, 


0 1  =  0° 
e2  =  45° 
e3  =  90° 
9i  =  180° 
05  =  315° 
e6  =  270° 


(23) 


Remark  5.1:  Note  that  the  placement  in  (23)  reflects 
the  placement  of  thrusters  about  the  USNA  yard  patrol 
experimental  platform  (see  Figure  4)  yet  the  barge 
dynamics  reflect  the  vessel  of  [17].  This  mistmatch  was 
allowed  since  accurate  drag  coefficients  have  not  been 
established  for  the  USNA  vessel. 

In  an  effort  to  limit  actuator  requirements,  the  fol¬ 
lowing  trapezoidal  velocity  desired  trajectory  [21]  was 
employed  for  point  to  point  movements/rotations 


Xd 


x0  +  f2 


Xf  +XQ  —  Vxtf 


0  <t<tb 

Vx  t  tb  <t<(tf  -  tb) 

x  f  ~  +  axt  ft  -  ft2  (tf  -tb)  <t<tf 


2 

axt% 


Fig.  2.  Simulated  Forces/Torques 


where  ax  and  Vx  represent  the  desired  barge  acceleration 
and  velocity,  respectively.  The  blending  time  tb  and 
the  arrival  time  tf  can  be  calculated  according  to  [21]. 
The  desired  position  and  rotation  trajectories  for  yd(t) 
and  ipd(t)  are  specified  in  an  analogous  manner.  For 
the  subsequent  simulation,  the  acceleration  and  velocity 
coefficients  were  selected  as 

ax  =  ay  =  0.1  (^),  =  0.0175  (^) 

Vx  =  Vy  =  \.Q  V^=0.1(^)  (24) 

Xf,Vf  =  100.0  (m),  ipf  =  0-0  (rad) 

For  the  command  vector  unet  of  (14)  and  the  force 
allocation  method  of  (20),  the  following  control  gains 
and  maximum  thruster  output  were  specified  as 

a  =  10.0,  kr  =  2.0,  Mmax  =  50.0  (N).  (25) 

Figure  2  illustrates  the  developed  forces/ torques  of  the 
swarm  configuration  from  the  optimization  allocation 
method  of  (20). 

VI.  Experimental  Results 

In  order  to  offer  experimental  validation  of  the  control 
strategies  presented  in  the  previous  sections,  a  small  scale 
experimental  vessel  was  built  (see  Figure  4).  The  hull  of 
the  “barge”  is  a  l/36th  scale  model  of  a  U.S.  Navy  Yard 
Patrol  craft.  It  measures  approximately  1.0  m  in  length 
and  0.3  m  wide.  Six  marine  bilge  pumps  were  attached 
to  the  vessel  to  act  as  our  swarm  of  tug  boats.  Each  bilge 
pump  is  powered  by  a  12.0  V,  H-bridge  power  amplifier, 
the  mechanical  design  of  the  pumps  is  such  that  they 
can  only  produce  thrust  in  one  direction.  The  maximum 
value  of  the  thrust  is  approximately  5  (N). 

The  control  strategy  is  executed  using  the  on-board 
Rabbit  3000  microprocessor.  Since  the  setup  is  intended 
to  be  utilized  inside  USNA’s  100.0  meter  tow  tank 
facility,  standard  global  positioning  sensors  could  not  be 
employed  to  provide  the  inertial  position  and  heading 
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(a)  (x,y)  Position  of  COM  of  Barge 
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Fig.  3.  (a)  ( x,y )  Position  of  Barge  COM,  (b)  Heading  ifj  of  Barge 


Fig.  4.  USNA  Experimental  Swarm  Vessel 


signals.  Instead,  a  standard  web  camera  was  mounted 
overhead  the  workspace  and  interfaced  with  MATLAB’s 
Image  Acquisition  Toolbox.  Two  high  intensity  LEDs 
mounted  on  top  of  the  vehicle  provide  active  color 
features  from  which  the  x,  y  position  and  heading  of 
the  vessel  in  the  image  plane  is  determined.  The  camera 
provides  a  640  x  480  pixel  image,  and  with  the  camera 
being  mounted  approximately  8.0  meters  overhead,  a  res¬ 
olution  of  approximately  0.02  m/pixel  is  obtained.  The 
position/orientation  measurements  were  transmitted  to 
the  ship-board  Rabbit  processor  via  pair  of  MaxStream 
Wireless  RS232  serial  modems.  The  required  velocity 
signals  of  P(t)  were  generated  through  a  backwards 
difference  scheme  coupled  with  a  low  pass  filter.  The 
numerical  values  for  the  mass  matrix  of  the  USNA  vessel 
differs  slightly  than  that  of  [17]  and  are  approximated 
as 

M  =  diag{  15.5  (kg),  15.5  (kg),  1.5  (kg  ■  to2)  } 

(26) 


x  (m) 


X  Desired 
X  Starting 
— ^  Measured 


(b)  Heading  Angle  of  Barge 


■  Desired 
1  Measured 


Fig.  5.  (a)  (x,y)  Barge  COM  Position  (b)  Heading  ijj  of  Barge. 


Since  exact  drag  coefficients  are  not  available  for  this 
vessel,  this  offers  a  perfect  scenario  under  which  to 
evaluate  the  adaptive  controller’s  ability  to  compensate 
for  unknown  drag  effects. 

Remark  6.1:  For  the  experimental  results,  we  would 
like  to  point  out  that  the  commutation  strategy  of  [2]  was 
utilized  to  allocate  forces/ torques  among  the  swarm  tugs 
in  lieu  of  the  optimization  strategy  of  (20).  In  addition, 
the  experimental  controller  operated  as  a  regulating 
controller  meaning  that  Pd(t)  =  Pd(t)  =  0). 

For  the  experiment,  the  barge  was  located  at  the 
following  initial  position/orientation  with  the  desired 
position/orientation  is  selected  as 

P(0)  =  [  4.64  (to),  0.93  (m),  205.40°  ]r 
Pd(t)  =  [  2.0  (m),  3.0  (to),  180.0°  ]T 

(27) 

The  control  gains  were  selected  as  follows 

a  =  0.3,  kr  =0.5,  (28) 

and  the  initial  values  for  the  parameter  estimate  vector 
0(f)  was  selected  as 

0(0)  =[  0.05,  1.5,  0.15  ]T  (29) 

The  (x,y)  position  of  the  center  of  mass  of  the  barge 
and  the  barge  heading  angle  V7  is  illustrated  in  Figure 
5  with  the  parameter  estimates  for  the  drag  effects  are 
shown  in  Figure  6. 

Remark  6.2:  It  should  be  noted  that  the  adaptive 
update  law  of  (15)  was  not  designed  for  identification 
of  the  drag  coefficients  D  =  diag {XU,YV,  Nr};  rather, 
the  dynamic  equation  of  (15)  was  specified  in  an  effort  to 
promote  the  position/orientation  tracking  error  objective 
by  ensuring  that  the  Lyapunov  function  V  (t)  is  negative 
semi-definite.  From  Figure  6,  one  can  observe  that 
the  parameter  estimates  approach  to  a  constant  value 
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Fig.  6.  Parameter  Estimates  for  Drag  Effects 


in  approximately  40.0  (sec)  which  coincides  with  the 
approach  of  the  disabled  vessel  to  its  desired  coordinates. 

VII.  Conclusions 

In  this  paper,  two  major  investigations  were  per¬ 
formed.  First,  an  optimization  based  force/ torque  alloca¬ 
tion  was  employed  and  compared  against  a  commutation 
based  force/torque  allocation  strategy.  This  optimization 
allocation  strategy  is  of  interest  due  to  the  commutation 
approach  being  susceptible  to  failure  in  heavy  sea  states. 
In  addition,  experimental  results  illustrating  the  ability 
of  the  adaptive  compensation  of  hydrodynamic  effects 
on  a  small  scale  experimental  test  stand  have  been 
illustrated.  Future  work  will  focus  on  the  inclusion  of 
force/ torque  mismatch  between  the  high  level  control 
unet  and  that  developed  by  the  collective  swarm  Bu 
within  the  stability  analysis.  Furthermore,  open-water 
experimental  verification  on  the  adjacent  Severn  River 
is  planned. 
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